Sound attenuation on the Bose-Einstein condensation of magnons in TlCuCl^ 
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We investigated experimentally and theoretically sound attenuation in the quantum spin system 
TICuGlz in magnetic fields at low temperatures. Near the point of Bose-Einstein condensation 
(BEC) of magnons a sharp peak in the sound attenuation is observed. The peak demonstrates a 
hysteresis as function of the magnetic field pointing to a first-order contribution to the transition. 
The sound damping has a Drude-like form arising as a result of hardcore magnon-magnon collisions. 
The strength of the coupling between lattice and magnons is estimated from the experimental data. 
The puzzling relationship between the transition temperature and the concentration of magnons is 
explained by their "relativistic" dispersion. 
PACS numbers: 75.30.Gw, 75.10.Jm, 78.30.-j 

Due to the richness and universality of phenomena that may be observed [1] quantum spin systems in magnetic 
fields have recently gained enormous attention. Similarities of quantum spin systems [2] and quantum fields allow to 
test experimentally and develop predictions of quantum field theories. The observation of plateaus in the high-field 
magnetization of one-dimensional quantum spin systems that follow the Oshikawa-Yamanaka- Affleck condition [3] is 
an example of this mutual interest. The effect of a magnetic field on a gapped spin system is especially interesting for 
a small gap, which can be reduced or even closed by experimentally available static fields. On the other hand sound 
attenuation experiments were very useful in the investigation of phase transitions with complex order parameters, e.g. 
superconductivity in the heavy fermion UPt3 [4] or the movement of flux lines in high-T c superconductors [5] . 

The XCUCI3 family of compounds, with X= Tl, K and NH4, realize especially interesting systems [2]. Magneti- 
zation plateaus are observed in NH4CUCI3, while TICUCI3 and KCUCI3 show field- induced critical phenomena. The 
underlying quantum magnetism is based on S — 1/2 spins of Cu 2+ ions arranged as planar dimers of CU2CI6. These 
dimers form infinite double chains parallel to the crystallographic a-axis of the monoclinic structure (space group 
• P2i/c with four formula units per unit cell) [6-8]. The ground state of TICUCI3 is a singlet with a spin gap A 7 K 
as a weakly anisotropic antiferromagnetic intra-dimer interaction of J ~60 K exists. The magnon dispersion of this 
■ compound was investigated theoretically in Ref. [9]. The three-dimensional (3D) character of the system leads to a 
relatively small gap since in this case it is more difficult to bind the magnons. In a magnetic field H the gap for 
S z = — 1 excitations is reduced as A — q/ibH where a g— factor varies between 2.06 and 2.23 depending on the H- 
' direction. Reaching a field of H g — A/g/iB, the gap is closed. 

A possible mechanism of magnetic-field induced spin transitions introduced by Giamarchi and Tsvelik [10] is the 
Bose condensation of the soft mode. For the explanation of field-induced magnetic ordering in TICUCI3 the Bose 
condensation of diluted magnons with S z = — 1 was proposed [11]. A Hartree-Fock description of such a transition 



[11] leads to an exponent $=1.5 relating concentration of magnons ng and the transition temperature Tg as rig ~ T 
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However, a larger value of $ cxp close to 2 has been determined by magnetization and specific heat experiments [11-14]. 

At the same time, the Bose condensation should be observable in experiments which are directly related to the 
distribution function of the magnons /s(fc), where k is the magnon momentum. Here we present, to our knowledge, the 
I first results of ultrasonic attenuation a(H) measurements on TICUCI3 in magnetic field and show that the experimental 
data suggests a gradual increase of fg(k) in the region of small k with the increase of H being consistent with the 
BEC of magnons at a critical field H c . As we will see below, ultrasonic attenuation experiments in this compound 
allow to trace the main features of the distribution function, investigate the fluctuation region where the attenuation 
anomalously increases, make a conclusion on the order of the transition, and determine the mechanism and strength 
of magnon-lattice coupling. 

Ultrasonic attenuation has been measured in single crystals of TICUCI3 prepared by the Bridgman method [12] 
using a pulse-echo method with longitudinal sound waves at a frequency of 5 MHz. The transducers were glued 
on freshly cleaved or on polished surfaces in the crystallographic (102) plane of the crystals using a liquid polymer 
(Thiokol, Dow Resin). The H- field has been applied parallel to the direction of sound propagation. The value of 
a(H) is determined from the experimental data as a(H) = 101g(//v/-T/v+i), where In and In+i are the intensities of 
Ath and Ath+1 pulses arriving the detector with the time interval 2t p , where t p is the pulse travelling time through 
the sample. 
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FIG. 1. Ultrasonic attenuation ot{H) with increasing (full circles) and decreasing magnetic field (open circles) and a fitting of 
the broad maximum by a Gaussian (dashed line). Inset (b): Transition regime after subtracting the Gaussian and a constant. 
Inset (c): Phase diagram with gH c /2 (thick line) from thermodynamic experiments [11]. The p-factor of S z = —1 excitations 
5 = 2.16. 



The experimental results shown in Fig. 1 can be summarized as the following observations: 

1) a weak iJ-dependence up to some upturn at H w 2.5 T, 

2) a sharp and asymmetric peak close to the critical field H c determined by thermodynamic experiments [11]. 
The sharp peak associated with the fluctuation region of the transition shows a pronounced hysteresis with increas- 
ing/decreasing field, which indicates a first-order contribution to the transition, and 

3) a decrease of a(H) for H > 9.7 T. 

If the comparably slow underlying increase is approximated by a Gaussian, the sharp peak near H c can be separated 
from the background. Inset (b) shows the resulting transition-induced attenuation at three different temperatures. 
The peaks have a gradual low- field onset of attenuation H° ns , a sharp maximum if™ ax and a high-field cut-off H^°. 
The two critical fields i?™ ax and H™ show a pronounced hysteresis. The respective attenuation at the critical fields 
systematically changes with T. For H° ns , that marks a crossover into a regime with strong fluctuations, the hysteresis 
is negligible. In a recent NMR investigation a hysteresis close to the phase boundary have also been observed and 
attributed to spin-phonon coupling [15]. Inset (c) shows the resulting phase diagram compared to magnetization 
measurements [11] (thick dashed line). 

The energy of a propagating sound pulse decays with time as: w(t) ~ exp[— (r + j m (H))t], and, therefore a(H) = 
201ge-[r + "f m {H)] t p . Here -fm(H) is the contribution of magnons coupled to the lattice, while an unknown background 
r arises due to all other effects like scattering by impurities, interfaces, etc. In contrast to F, the -ff-dependence of the 
damping ~f m (H) — 7 m (0) which corresponds to the field-induced changes in and sound attenuation mechanism 

can be determined accurately from Fig. 1. By comparison with theoretical calculations also 7 m (0) will be obtained 
from the experiment. 

Cottam [16] investigated spin-phonon coupling in an antiferromagnet in the case of a low density of magnons where 
magnon-magnon interactions can be neglected. Below we investigated another case where the damping is related 
to collisions between magnons. The model for sound attenuation is the following: In the field of the longitudinal 
sound wave with displacements of ions u(r), the energy of the magnon changes due to magnon-lattice coupling as 
<5e(k) = CfeVu(r), where Ck is the deformation potential. The spatial dependence of the displacement in a plane wave 
with wavevector q and frequency u>: u(r) = Uo cxp[z(qr — ut)] leads to an effective time-dependent force F = — Vfe(k) 
with F = CkUoq 2 driving changes in the velocity of magnons. Since one cannot create two S z = — 1 states on a Cu-Cu 
bond, the magnons collide like hard-core bosons. Inelastic collisions between the magnons driven by the displacements 
of the lattice in the wave lead to a dissipation of the pulse energy. The dissipation rate per one magnon is equal to 
the average of the product Fv per period, where v is the magnon velocity. The sound dissipation rate in this case is 
the dissipation rate of an external harmonic field in a medium, generally described by a Drude-like mechanism of the 
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form: 
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where nik is the magnon effective mass and i>k.k' is the relative velocity of the particles with momenta k and k' that 
determines their collision rate. The constant A is of the order of one and depends on details of magnon-magnon 
collisions, which are beyond our consideration. The mean time between the collisions (r) ~ {iv^ 1 ) with the magnon 
free path £ being proportional to l/na, where a is the magnon-magnon scattering cross-section, and n is the magnon 
concentration, is, therefore sensitive to the distribution function. The mean value (v^ 1 ) and the damping increase 
simultaneously. Since cj(t) is very small (<~ 0.01) at the experimental u) ~ 3 • 10 7 s _1 , it will be neglected in the 
denominator of Eq.(l) for the rest of the discussion. 

To include the specific spin excitation spectrum, we start with the " relativistic" form of the magnon dispersion 
e(fc) = y/ A 2 + J 2 fc 2 /4 at small k which is important if Tb and A are of the same order of magnitude as in TICUCI3. 
From here on the lattice constant I = 1 if otherwise is not stated. Since thermal k <C 1 at Tb <C J this form of energy 
is sufficient for an understanding of the Bose condensation in the experimental temperature range. The velocity of a 
magnon and its inverse mass are (we put U = 1): 
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respectively. The mass at small k, m — 4A/J 2 , determines Tb at low concentration of magnons where Tb <C A. In 
the mean-field approximation which will be used below, the distribution function is Jb{s) = 1/[ CX P((^ — l^cs)/T) — 1], 
where e{k) = e{k) — A, and ^ g' = g^siH — H g ) — 2V lam n. n is the dimensionless concentration of the magnons 
determined by the magnetization MasS= M/H, and V mm is the magnon-magnon interaction parameter. The 
transition temperature Tb obtained with the condition fi e g = 0, that is 
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(3) 



is, therefore, Vmm-independent. Fig. 2a presents dimensionless Hb{Tb), and Fig. 2b shows the ratio 
(Tlb{Tb) /nsi^)) /Tg obtained from Eq.(3) for different A and J = 7A. As one can see, at the chosen A = 6.5 K, 
the ratio is nearly a constant in the experimental range of temperatures, giving an explanation for the observation of 
(f) cxp rts 2 [11-14]. Since near the condensation point the thermal wavevector of the magnon k 2 ~ mTs, leading to a 
temperature dependence of the mean value (e(k) — A) ~ T, the "exponent" depends on Tb/ A only. We note that 4> exp 
is indeed a fitting parameter, which can depend on the measured quantity. We also mention that the parameter which 
determines the applicability of the Drude approach ki and can be estimated as <~ n~ 2//3 , where k is the characteristic 
magnon momentum, is always larger than ten thus justifying the validity of Eq.(l) except the fluctuation region. 
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FIG. 2. (a) Dimensionless concentration of magnons ub(Tb) corresponding to the BEC temperature Tb- Squares present 



experimental data 
corresponding to <f> 



[11]. (b) [n B (T B )/n B (IK)] T~' 
= 3/2. 



With increasing A the curves are closer to T B (the dashed line), 



Since near the transition the Bose function behaves as 



fa{e) 
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(4) 



due to the singularity at ji e e — > 0, the average diverges as ln(TB/|/i c ff|) leading, as it will be shown below, to 

an increase in the damping rate. 

The sound-induced change in the magnon energy consists of two terms, fe(k) = 5ea0*) + fej(k) arising due to a 
modulation in the gap A = A + CaVu and the exchange J = J + CjVu, respectively. The resulting deformation 
potential is: 
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Since the J-originated term in Ck vanishes at small k, the main contribution to the increase of damping near the 
transition comes from Ca- The calculated behavior of the damping at T=3.5 K for H < H c is shown in Fig. 3 for two 
VInm and two models of phonon-magnon coupling. As a result of the increase of /s(fc) at small k near the transition 
point, the damping due to the C*A-term increases and shows a 7 m (lZ) behavior which is in full agreement with the 
experimental data in Fig. 1. At the same time, due to a vanishing coupling for k — > in the second contribution 
in Eq.(5), the Cj— dependent term leads to a slight decrease in the damping, in contrast to Fig.l. Moreover, near 
the transition point, where k 2 ~ TbA/J 2 , the ratio of the Cj— and Ca— originated terms Eq.(5) is of the order of 
CjTb/JCa, thus containing a small factor of Tb/J that decreases the role of the Cj— originated term. 
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FIG. 3. Calculated fag (a) and damping (b,c) for different spin-lattice couplings as function of field H at T — 3.5 K. For 
H < 2.5 T the behavior of damping is almost flat, in agreement with the experiment. A = 6.5 K, H g = 4.5 T. The results 
depend on Vinm which determines the H c . 

Our approach can qualitatively explain the decrease in the damping rate when the field exceeds some value cor- 
responding to the broad maximum at H = 9.7 T in Fig. 1. The above consideration is based on the low-energy 
singularity of the /s(e) increasing as 1/e at H — H c . At H > H c the number of quasiparticles n rapidly increases 
with H [11]. The spectrum of excitations in this case is E(k) — v c k, with v c ~ n 1 / 2 being the sound velocity in the 
condensate. With the increase of H, and, correspondingly n, the excitation energy increases, and, therefore, the role of 
thermal excitations decreases. Thus, the increase in H drives the condensate effectively towards the zero-temperature 
behavior. The system demonstrates a zero-T behavior when n becomes much larger than ns corresponding to the 
BEC at given T. At zero T the distribution of the out-of condensate particles which contribute to the damping at 
small momenta is Ne(E — > 0) = mv^/2E(k) [21]. Since E(k) is linear in k, the singularity in Eq.(4) at // e ff = 
becomes weaker, and the damping decreases. To estimate the field at which the behavior of the condensate is close to 
the T = limit we note that at T > thermal excitations and interaction V^ lm produce out-of condensate particles. 
The concentration of the magnons at which the system demonstrates a crossover to the T — behavior is determined 
by the condition that the concentration of the particles out of the condensate arising due the interaction (~ n\/na 3 ) 
is larger than the concentration of thermally excited particles proportional to T 3 / 2 , where a = V mm ? 3 m/47r is the 
magnon-magnon scattering length. In TICUCI3 a and I are close to each other since the hard-core boson has a spatial 
size of the Cu-Cu distance. As a result, the critical concentration which determines the crossover to zero-temperature 

1 /3 

behavior and, in turn, the broad maximum position is n cr ~ TiB\l I 'an B ) ~ lOng. On the other hand from the data 
on n(H)— dependence of Ref. [11] one expects that at T = 3.5 K, n(H = 10T) is an order of magnitude larger than 
nB{E[ c ), in a good qualitative agreement with our estimate. 
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Following Eqs.(l) and (5) one can estimate Ca- The energy density in the sound wave with ui = cq (c is the sound 
velocity) is: w ~ kuqlu 2 /c 2 , where k is the elastic constant. The dissipation rate estimated from Eq.(l) is: 

dw F 2 1 2 ^ 4 2 1 

-77 rr-^A— «o 7 r • (6) 

ar m (T (u) c 1 ma [vj 

Therefore, the relaxation rate (dw/dtyw^ 1 <~ w 2 . Let us consider as an example 7 m (0) that with the estimates given 
above can be written as: 

, n . dwl C 2 A cu 2 l J 

where we took into account that far from the transition point (l/v) ~ 1/ (v) with (v) ~ \JT jm. From the theoretical 
7 m (_ff)— dependence (Fig. 3), we conclude that j m (H g )/ 7 m (0) <~ 2. From the experimental data in Fig. 1 we obtain 
j m (H g )— 7 m (0) ~ 0.1 dB, and, therefore, 7 m (0) corresponds to a damping of the order of 0.1 dB. For the thickness of 
the sample of 2 mm and the sound velocity c = 3 • 10 5 cm/s, the pulse travelling time t p is 6.6- 10 -7 s, and, therefore, 
7m (0) ~ 0.2 • 10 4 s _1 . The density p — 5 g/cm 3 , ui ~ 3 • 10 7 s -1 , and / = fi 1 / 3 w 5 A, where $7 is the unit cell volume 
per chemical formula yield Ca of the order of 200 K. 

We presented ultrasonic attenuation experiments as well as a theoretical analysis on the field-induced changes in 
the spin dimcr system TICUCI3. These experiments are in agreement with the magnon Bose-Einstein condensation 
and summarized as a sharp peak in the attenuation a(H) in the close vicinity of the transition and a second broader 
maximum at larger fields, deep in the long-range ordered phase. The peak in a{H) shows a hysteresis crossing the 
phase boundary with increasing/decreasing field. This effect is attributed to a first-order contribution to the phase 
transition. The exponent <f> w 2 in the relation hb ~ Tg and the iJ-dependence of sound attenuation originate 
from the 'relativistic' dispersion of magnons. The spin-lattice coupling is due to a phonon-induced modulation of the 
spin gap with the coupling constant Ca ~ 200 K. The recent observation [22,23] of triplet crystallization in the 2D 
compound SrCu2(BC>3)2 in large magnetic fields with similar aspects of spin-phonon coupling demonstrate that these 
phenomena are of a general nature. 

We acknowledge important discussions with D. Lcnz, M. Takigawa, B. Wolf and B. Luthi, and financial support by 
DFG/SPP 1073 and the Austrian Science Fund project P-15520. Experiments have been performed at the High-Field 
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